Density-dependent effective nucleon-nucleon interaction 
from chiral three-nucleon forces* 

J. W. Holt, N. Kaiser, and W. Weise 
Physik Department, Technische Universitat Miinchen, D-85747 Garching, Germany 

Abstract 

We derive density-dependent corrections to the in-medium nucleon-nucleon interaction from the 
leading-order chiral three-nucleon force. To this order there are six distinct one-loop diagrams con- 
tributing to the in-medium nucleon-nucleon scattering T-matrix. Analytic expressions are presented 
for each of these in both isospin-symmetric nuclear matter as well as nuclear matter with a small 
isospin asymmetry. The results are combined with the low-momentum nucleon-nucleon potential 
Vj ow _k to obtain an effective density-dependent interaction suitable for nuclear structure calculations. 
The in-medium interaction is decomposed into partial waves up to orbital angular momentum L = 2. 
Our results should be particularly useful in calculations where an exact treatment of the chiral three- 
nucleon force would otherwise be computationally prohibitive. 



* Work supported in part by BMBF, GSI and by the DFG cluster of excellence: Origin and Structure of the 
Universe. 



I. INTRODUCTION 



Three-body forces play an essential role in describing the detailed properties of nuclear 
many-body systems. Compelling evidence for this arises from the inability of non-relativistic 
nucleon-nucleon (NN) interactions alone to accurately reproduce (i) the saturation properties of 
nuclear matter, (ii) the binding energies and spectra of light nuclei, and (iii) nucleon-deuteron 
scattering differential cross sections at intermediate energies. These systems are simple enough 
that the uncertainties associated with the solution of the corresponding many-body equation are 
well controlled, and therefore deviations from experimental data are interpreted as deficiencies 
in the underlying nuclear force model. 



In the case of nuc 
hole-line expansion 



e^r matter, a variety of methods such as the Brueckner-Bethe-Goldstone 
variational methods 3|,|4j, and the summation of particle-particle hole- 
hole ring diagrams [5|, |6|] all produce similar results for the energy per particle as a function of 
the density. However, these calculations result in a series of saturation points that lie at either 
too high a density or too low a binding energy compared to the empirical values of po — 0.16 
fm -3 and E/A = —16 MeV. Numerous solutions to this problem have been proposed, such 



as includin g re lativistic effects 



masses 



11 



three-nucleon forces 



101 ]. or medium-modified meson 



121 ] . The analogous problem in light nuclei is that calculated binding energies come 



out ~ 10 — 15% smaller than experimental values when two-body interactions alone are used 



14 



15] . By including various models of the three-nucleon force, several groups have shown 



em can be remedied but the low-lying spectra of 



16 



171 ] . Finally, the differential cross section for 



that not only the nuclear binding energy prob 
light nuclei can also be greatly improved 
elastic nucleon-deuteron scattering at intermediate energies {Effi = 100 — 300 MeV) exhibits 
a minimum at backward angles that is larger by nearly 30% compared to Faddeev calculations 
using NN interactions alone ljj, 3- This discrepancy is well accounted for when effects from 



three-body forces 



20] or (equivalently) explicit A isobar degrees of freedom lg, |2l| are included. 



The preceding examples illustrate the limitations in obtaining precise agreement with var- 
ious nuclear observables when using only two-body interactions that are fit to free-space NN 
scattering data. Although there has been much effort devoted to the construction of realistic 
models of the three-nucleon force, it is technically very challenging to include them in large-scale 
calculations of medium and heavy nuclei. An alternative and simpler approach is to employ in- 
stead density-dependent two-body interactions that reflect the underlying three-nucleon forces. 
Therefore, the purpose of the present work is to construct a density-dependent in-medium 
NN interaction that is generated at one-loop order by the genuine three-nucleon force. We 
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have found in ref. |22J that such density-dependent corrections to the NN interaction strongly 
suppress the 14 C ground state to 14 N ground state Gamow- Teller transition and in this way 
helps to explain the anomalously long lifetime of 14 C. We work within the framework of chiral 
effective field theory and take into account the leading-order three-nucleon forces arising in 
that systematic approach. We consider first the density-dependent terms in isospin symmetric 
nuclear matter and then continue with the leading (linear) effects due to an isospin asymmetry, 
which are likely to be relevant for heavy nuclei. The paper is organized as follows. In section 
|II]we briefly describe the chiral effective field theory approach to nuclear interactions, together 
with the explicit form of the leading-order three-nucleon force. In section II III we evaluate the 
one-loop corrections to NN scattering in the nuclear medium generated by this chiral three- 
nucleon force. In section [IV] we present numerical results for the in-medium NN interaction in 
partial waves up to orbital angular momentum L = 2. The paper ends with a summary and 
conclusions. 

II. NUCLEAR INTERACTIONS FROM EFFECTIVE FIELD THEORY 
A. Chiral effective interactions 

Chiral effective field theory provides a useful tool for describing a wide range of low-energy 
hadronic phenomena within a single consistent framework. The theory exploits the natural 
separation of energy scales that results from the spontaneous (and explicit) breaking of chiral 
symmetry in QCD, which gives rise to light pseudoscalar Goldstone bosons (the pion triplet 
in the two flavor case). These, together with nucleons, comprise the low energy degrees of 
freedom of the theory. Long-range effects from one- and two-pion exchange between nucleons 
are treated explicitly, while the short- distance dynamics due to heavy mesons and baryon 
resonances (with the possible exception of the A isobar) are integrated out and their effects 
are encoded in nucleon contact terms. Presently the computation of the chiral NN potential 
has been carried out to next-to-next-to-next-to-leading order (N 3 LO) in the small momentum 



expansion 



23 



24 



25 



261 ] . At this order it is possible to achieve an agreement with empirical NN 



scattering phase shifts that is comparable to previous high-precision NN potentials [27], 1281. |29| . 
By adjusting the 29 parameters (mostly low-energy constants) that occur at this order, the 
1999 database for np and pp elastic scattering up to = 290 MeV can be fit with a x 2 /dof of 
1.1 for np scattering and 1.5 for pp scattering. Furthermore, the experimental deuteron binding 
energy, charge radius, and electric quadrupole moment are very well reproduced by the chiral 
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N 3 LO potential 



1 



When applied to two- and few-nucleon problems, these chiral potentials 



are regulated by exponential functions 



23, 



30t ] with cutoffs ranging from 500 to 700 MeV in 



order to eliminate high-momentum components. In addition, in-medium chiral perturbation 
theory which emphasizes the role of explicit two-pion-exchange dynamics in nuclear matter has 
been developed in ref. [3lfl . 



Three-nucleon forces arise first at third-order in the chiral power counting |32j. Three com- 
ponents of different range, shown by diagrams (a),(b), and (c) in Fig. [TJ, occur at this order and 
have the following analytic structure 
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where cfi = pi' — pi is the difference between the final and initial momentum of nucleon i and 



Ffh = ^ Hcml + 2c 3 % • qj) + c^rla k ■ (gj 



x 



(4) 



Together with = 1.29 and f n = 92.4 MeV, the parameters of namely c\ 



0.81 GeV 1 , C3 = — 3.2 GeV 1 , and C4 = 5.4 GeV 1 , are well determined from fits to low- 



energy NN phase shifts and mixing angles 



231 ] . Restricting the two large coefficients c 3)4 to 



their dominant A(1232)-resonance contributions, one reproduces the celebrated three-nucleon 



force of Fujita and Miyazawa [33j. The medium-range (V^ ) and short-range (V^ } ) com- 
ponents are proportional to two new low-energy constants and ce, respectively. These 
constants can be fixed by fitting properties of few-nucleon systems, such as the triton binding 



r(cth 



energy together with the nd doublet scattering length 



341 ]. the 4 He binding energy [35j, or 



binding energies and spectra of light nuclei [17J. At next order (in the chiral expansion) there 



32], but 



are many additional one-loop diagrams contributing to the chiral three-nucleon force 
no new low-energy constants appear. The explicit construction of all these terms is currently 
in progress. 



B. Low-momentum nuclear interaction 



Recently there has been much interest in understanding the scale-dependence of NN inter- 
actions from the point of view of the renormalization group 



36 



371 ] . In most applications of 
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fa) 



fb) 



(c) 



FIG. 1: The leading-order contributions to the chiral three-nucleon force: (a) the long-range 2tt 
exchange force V^j ■ , (b) the medium-range l7r-exchange force and (c) the short-range contact 



interaction K 



(ct) 



37V • 



chiral effective field theory to few-nucleon systems, exponential regulator functions with cutoffs 
between 500 and 700 MeV are used to eliminate the high-momentum components from the 
NN interaction. In fact, all realistic NN potentials are fit to pp and pn scattering data below 



a laboratory energy of 350 MeV and therefore are constrained experimental 
momentum of p max — 400 MeV = 2.1 fin -1 . Bogner and collaborators 36 



y 



only up to a 



37l | have shown 



how to evolve any bare NN interaction down to a low-momentum scale via renormalization 
group techniques. Such low-momentum interactions Vi ow -k are phase-shift equivalent to the 
underlying bare interaction at energies below the specified cutoff scale Ai ow _k. Moreover, when 
the decimation scale is reduced to p max , all realistic NN interactions merge to a nearly universal 
potential, thereby removing the model dependence related to the high momentum components. 
The method for constructing such low-momentum interactions is as follows. 

One begins with the half-on-shell T-matrix for free NN scattering in a given partial wave, 



Tip ,p) = V NN {p ,p) + -VI 5 5 q dq, 

K Jo 



p2 _ g2 

and introduces a low-momentum half-on-shell T-matrix at the scale A 



low— k 



T low _ k (p,p) = Viow-k (p' ,p) + -V 



Vio W -k(p', g)T low _ k (g,p) 2 



p2 _ g2 



(5) 



(6) 



JO 

where V denotes the principal value. Since one has to preserve low-energy observables, e.g. 
scattering phase shifts defined through tan5(p) = —pT(p,p), one requires that 



Ti ow _ k (p',p) = T(p',p) , for p',p < Aiow-k- 



(7) 



Eqs. [jp and (JTj) together define the low-momentum NN interaction Vi ow _k (for a review see 
ref. [37j]). Under the scale decimation procedure, all high-precision NN potentials flow, as 



Aiow-k — ► 2.1 fm _1 , to a nearly unique low-momentum potential Vj ow -k- In the present study we 
employ the (bare) chiral N 3 LO interaction and evolve it down to a resolution scale Ai ow _k = 2.1 
fm _1 . 



As demonstrated in ref. [10] the addition of three-nucleon forces is essential for obtaining 
reasonable saturation properties of nuclear matter when using the universal low-momentum NN 
potential V[ ow _k in Hartree-Fock calculations. The construction of decimated low-momentum 
three-body forces consistent with the two-body decimation for Vj ow -k remains a challenge. The 
difficulty has so far been addressed pragmatically by exploiting that low-energy nuclear observ- 
ables must be scale- independent. Following this reasoning the low-energy constants associated 
with the one-pion exchange component c,D(Ai ow _k) and the short-range contact term ce(Ai ow _k) 
of the chiral three-nucleon interaction have been fitted to experimental binding energies of three- 
and four-nucleon systems ( 3 H, 3 He, and 4 He) at a given (variable) low-momentum scale Ai ow _k- 
Since we investigate the density-dependent two-nucleon interaction at the scale Ai ow _k = 2.1 
fm _1 , we use the corresponding values for cp and c# obtained in ref. 



c D (2.1 far 1 ) = -2.062 , c E (2.1 far 1 ) = -0.625 , (8) 

together with A x = 0.7 GeV. In addition to the low-momentum decimation techniques de- 
scribed above, there exist several other means to construct effective nucleon-nucleon interac- 
tions suitable for perturbative many-body calculations, such as the similarity renormalization 



group transformations [39( and the unitary correlation operator method (UCOM) [40j, l4lj |. 
The resulting effective two-body interactions are all quantitatively similar, and the inclusion of 
three-nucleon forces within these different frameworks is currently in progress. 



III. IN-MEDIUM NUCLEON-NUCLEON INTERACTION 

A. Density-dependent terms in isospin-symmetric nuclear matter 

In this section we derive from the leading-order chiral three-nucleon interaction, eqs. ([Ti- 
ll]), an effective density-dependent in-medium NN interaction. We are considering the on- 
shell scattering of two nucleons in isospin-symmetric (spin-saturated) nuclear matter of density 
p = 2k 3 j-/37T 2 in the center-of-mass frame, Ni(p) + A^— p) — * N\(p + q) + A^— p — q). The 
kinematics is such that the total momentum of the two-nucleon system is zero in the nuclear 
matter rest frame before and after the scattering. The magnitude of the in- and out-going 
nucleon momenta is p = \p \ = \p + q\, and q = \q\ is the magnitude of the momentum transfer. 
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T 



(3) 



FIG. 2: In-medium NN interaction generated by the two-pion exchange component (~ ci,3 4) of the 
chiral three-nucleon interaction. The short double-line symbolizes the filled Fermi sea of nucleons, i.e. 
the medium insertion — 2-7r<5(£;o) &(kf — \k |) in the in-medium nucleon propagator. Reflected diagrams 
are not shown. 

The on-shell interaction in momentum-space has the following (general) form 

V(p, q) = Vc + Ti ■ r 2 W c + [V s + f\ ■ f 2 W s ] <?i ■ <?2 + [Vt + n ■ t 2 W t ] fi\ ■ qo 2 ■ q 
+ [Vso + ti • f 2 W so ] + ^2) ■ (qxp) 

+ [Vq +n-r 2 W Q ] a x - (qxp)a 2 - {qXp). (9) 

The subscripts refer to the central (C), spin-spin (S), tensor (T), spin-orbit (SO), and quadratic 
spin orbit (Q) components of the NN interaction, each of which occurs in an isoscalar (V) and 
an isovector (W) version. For the purpose of comparison with the density-dependent terms 
that follow, we reproduce the expression for the (bare) l7r-exchange: 

(lvr) _ 9^Mh 



V, 



NN 



167T/2 



n -> - 0-1 ■ q c 2 • q 

T\ ■ T 2 



mi 



(10) 



Here a\^ 2 and 7^2 are the usual spin and isospin operators of the two nucleons. Note that 
we have included an additional factor of Mn/4tt in Vnn in order to be consistent with the 
conventions commonly chosen for V[ ow -k- 

We start with those contributions to the in-medium NN-interaction V^ d that are generated 
by the 27r-exchange component of the chiral three-nucleon force. The three different topologies 
for non-vanishing one-loop diagrams are shown in Fig. [2j The short double-line on a nucleon 
propagator symbolizes the filled Fermi sea of nucleons, which introduces the medium insertion 
— 27r5(k ) 8(kf — \k\) in the in-medium nucleon propagator. In effect, the medium insertion 
sums up hole propagation and the absence of particle propagation below the Fermi surface 
\k I < kf. The left diagram in Fig. [2] represents a l7r-exchange with a Pauli blocked in-medium 
pion self-energy and the corresponding contribution to V^ d reads: 



mcd,l 9A M NP _ _ &l ■ ■ q , n 2 , 2^ 



NN 



T\ ■ T 2 



ml + q 2 ) 2 



(2dm n + c z q A 






FIG. 3: In-medium NN interaction generated by the one-pion exchange (~ Co) and short-range 
component (~ ce) of the chiral three-nucleon interaction. 

Since c\$ < 0, this term corresponds to an enhancement of the bare l7r-exchange. It can be 
interpreted in terms of the reduced pion decay constant, f* 2 s — f% + 2c 3 p, that replaces f 2 in 
the denominator of eq. (fTUl) and which is associated with the space components of the axial 
current in the nuclear medium. The second diagram in Fig. [2] includes vertex corrections to the 
l7r-exchange caused by Pauli blocking in the nuclear medium. The corresponding contribution 
to the in-medium NN-interaction has the form: 



v 



med,2 
NN 



in _ a 1 -qa 2 -q i 2 
rri • r 2 5 ; <J - Aam^ 



T (p) + T 1 (p) 



32tt 3 /^ ml + q 2 

-(c 3 + c 4 ) \q 2 (^o(p) + 2r 1 ( P ) + r 3 (p) ) + 4r 2 ( P ) 



+ 4c 4 



2& * n 

Y-mJro(p) 



.(12) 



Here, we have introduced the ^/-dependent functions r.,(p) which result from Fermi sphere 
integrals, f^ <kf d 3 k, over a static pion-propagator [ml + (k + p) 2 ]^ 1 : 



r oG°) =kf — m jr 

r iW = +p 2 )-r (p) 
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arctan — h arctan — 
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^(p) = ? + g(*J 



p 2 )Tv(p) + -(ml + k 2 -p 2 )T x (p) 



r 8 (p) 



3p 2 



m 2 + fc 2 + P 2 

V 



r (p) 



fc 2 + 3p 2 



2p 2 



ri(p). 
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(14) 

(15) 
(16) 



By analyzing the momentum and density dependent factor in eq. (fT2|) relative to V^ } , one finds 
that this contribution corresponds to a reduction of the l7r-exchange in the nuclear medium. 
Approximately, this feature can be interpreted in terms of a reduced nucleon axial- vector con- 
stant g* A . 

The right diagram in Fig. [2] represents Pauli blocking effects on chiral 2-7r-exchange. Evaluat- 
ing it together with the reflected diagram one finds the following contribution to the in-medium 



NN-interaction: 



V, 



med,3 
NN 



9aM. 



N 



64vr3/ 4 



12 Cl ml 



2T (p)~(2m 2 n + q 2 )G (p,q) 



-c 3 



8kj - 12(2/7^ + q 2 )T (p) - Qq 2 T l (p) + 3(2m 2 7r + q 2 ) 2 G (p, q) 



+4c 4 f! ■ r 2 (<Ti ■ a 2 q 2 - #i ■ ga a • q)G 2 (p, q) 

-(3c 3 + c 4 fi • f 2 ) zfa + a 2 ) • (g x p) [2r (p) + 21^) - (2ml + tf) 
x (G (p, q) + 2Gi(p, g))] - 12 Cl m2 + a 2 ) ■(gxp) [g (p, g) + 2Gi(p, g) j 

+4c 4 f 1 -f 2 o i i-(gxp)a : 2-(gxp) G (p,g)+4Gi(p,g)+4G 3 (p,g) J. (17) 

One observes that in comparison to the analogous 27r-exchange interaction in vacuum (see 
section 4.2 in ref. (42J) the Pauli blocking in the nuclear medium has generated additional 
spin-orbit terms, + a 2 ) ■ (q x p), and quadratic spin-orbit terms, a\ ■ (qx p) a 2 ■ (qx p), 
written in the last three lines of eq. ([IT]) . The density dependent spin-orbit terms (scaling with 
03,4) in the in-medium NN-interaction V^ d ' 3 demonstrate clearly and explicitly the mechanism 



of three-body induced spin-orbit forces proposed long ago by Fujita and Miyazawa [33(. The 
functions Gj (p, q) appearing in eq. ( ITTl) result from Fermi sphere integrals over the product of 
two different pion-propagators. Performing the angular integrations analytically one arrives at: 



Go,*,**(p, g) 



dk 



{k,k 3 ,k 5 } , qk+ JA(p) + q 2 k 2 
In 



18) 



QJo y / A(p) + qW y/A(p) 

with the abbreviation A(p) = [ml + (k +p) 2 } [m 2 + (k — p) 2 ] . The other functions with j = 1,2,3 
are obtained by solving a system of linear equations: 

r„(p) - {ml + p 2 )G (p, q) - G*(p, q) 



Gi(p,q) 
Gu(p,q) 



4p 2 — q 2 
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G 2 (p,q) = (ml + p 2 )G 1 (p,q) + G*(p,q) + G u (p,q), 
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(20) 
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In this chain of equations the functions indexed with an asterisk play only an auxiliary role 
for the construction of £71,2,3 (p, g). We note that all functions Gj(p,q) are non-singular at 
q = 2p (corresponding to scattering in backward direction). For notational simplicity, the 
^/-dependence of Tj(jp) and Gj(p, q) has been suppressed. 

Next, we come to the l7r-exchange component of the chiral three-nucleon interaction pro- 
portional to the parameter cd/A x , where cp ~ —2 for a scale of A x = 0.7 GeV [38j. The filled 
black square in the first and second diagram of Fig. [3] symbolizes the corresponding two-nucleon 
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one-pion contact interaction. By closing a nucleon line at the contact vertex, one obtains a 
vertex correction (linear in the density p) to the l7r-exchange: 



V, 



med,4 
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g A M N c D p _ _ ai-qa 2 -q 
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(23) 



32tt/*A x * - rrii + q' 
Since cd is negative, this term reduces again the bare l7r-exchange, roughly by about 16% 
at normal nuclear matter density po — 0.16 fm -3 . The second diagram in Fig. [3] includes 
Pauli blocked (pionic) vertex corrections to the short-range NN interaction. The corresponding 
contribution to the density dependent in-medium NN interaction reads: 



V 'Hied, 5 
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64tt3/4A x 



—,^1 -(qxp)a 2 -(qx p) 
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(24) 



where we have used an identity for 3\ • p a 2 • p + a\ • (p + q) a 2 • (p + q) = [...]. The ellipses 
stands for the combination of spin operators written in the square bracket of eq. (T2"4"j) . 

Finally, there is the short-range component of the chiral 3N interaction, represented by a 
three-nucleon contact- vertex proportional to c^/A x . By closing one nucleon line (see right 
diagram in Fig. 2) one obtains the following contribution to the in-medium NN-interaction: 

med,6 

3M N c E p 



V 



NN 



(25) 



which simply grows 



nucleon momenta 44]. 



inearly in density p = 2k 3 j-/3ir 2 and is independent of spin, isospin and 



The above expressions have been derived for on-shell scattering, which greatly simplifies 
the calculation of the density-dependent corrections to the NN interaction. As discussed 
in ref. 22], a suitable choice for extrapolating these expressions off-shell is to make the 
substitution p 2 — > |(p 2 + p' 2 ). Additionally we include a regulator function of the form 
exp [— (p/Ai ow _k) 4 — (p',Miow-k) 4 ], where A[ ow _k = 2.1 fm -1 is chosen so that the two-body 
and three-body force contributions are decimated down to the same scale. This is consistent 
with the approach taken in ref. 38]. 



B. Correction terms in isospin asymmetric nuclear matter 

Now we consider the additional modifications W^jy' 1 to the six in-medium NN scattering 
T-matrices, V^ d '\ resulting from a small isospin asymmetry. The total nucleon density made 
up by protons and neutrons is p = p n + p p , and the Fermi momentum kf is given as before 
by the relation p = 2/c 3 c /37r 2 . As a measure of the isospin asymmetry, we define the relative 
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neutron excess 5 np = (p„ — p p )/ p. For a heavy nucleus such as 208 Pb, the relative neutron excess 
is of the order 5 np ~ 0.2. Re-evaluating the six diagrams in Figs. [2] and [3] with the substitution 



0(k f -\k\) 
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' np 3 dk 



(26) 



we obtain the following expressions for the corrections (linear in 5 np ) to the density-dependent 
NN interaction. 
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4cim„ 



2r (;p)-(2m 2 r + g 2 )G'o(p,g) 



-c 3 



3^3 

l _£ _ 4(2m 2 + g 2 )ro(p) _ 2 g2 ri(p) + (2m 2 + g 2 ) 2 G (p, g) 



+4c 4 (<?i -a 2 q 2 qa 2 - q)G 2 (p, q) 

-(c 3 + c 4 ) ifa + <7 2 ) • (g x p) ^2T (p) + 2T 1 (p) - (2m 2 + q 2 ) 

x(G (p,q) + 2G 1 (p,q))] - 4cim 2 i(Bx + B 2 ) ■ (qX p ) \g (p, q) + 2Gi(p, q) 



+4c 4 ai ■ (g x p)B 2 - (?x p) G (p, g) + 4Gi(p, g) + 4G 3 (p, g) 



(29) 



med,4 

g A M N c D (p p Pn 



NN 



7774 + g 



(30) 



med,5 QaMnCd^up , 3 3% , 1 - - 1 , . 



384tt 3 /4A x 

f 2p 2 \ 2 
x M - — J - — B x ■ (g x p) o 2 ■ (g x p) 



B\ • (To( 2p 2 - y 



r (p) + 2r 1 (p) + r 3 (p) 



a x -qa 2 -q 



and (31) 



W 



med,6 

3M N c E (p p - p n ) 



NN 



16tt/4A, 



+ rf ) . 



(32) 



11 



For completeness we also list the derivatives with respect to kf of the loop functions Tj(p) and 
Gj(p,q) encountered in the previous section: 

dT (p) 



dT 2 (p) 
dk f 



dk f 

dTijp) 
dk f 

k 2 

^(ml + kj+p 2 ) 

dT 3 (p) 
dk f 
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21 

p 2 2p 2 dkf ' 

- — 2 K + (*/ + p) 2 } H + (k f - P?] ^ 



dk 



f 



1 

p 2 



,dr (p) dv 2 {p) 



*' dkf 



-3- 



dk f 



_d_ 

dk t 



2{kf,k 2 ,k 3 f} i qkf + ^A(p) + q 2 k 2 



(33) 
(34) 
(35) 
(36) 

(37) 



qy/A(p) + q 2 k) y/A(p) 
with the abbreviation A(p) = [m 2 + (kf + p) 2 ] [ml + (kf — p) 2 ]. Aside from the pion self- 
energy term, which receives no modification linear in the isospin asymmetry, each term W^ffi 1 
is proportional to the factor (rf + rf ) 5 np . Therefore, the proton-proton (pp) and neutron- 
neutron (nn) interactions in isospin-asymmetric nuclear matter receive corrections of opposite 
sign, while the neutron-proton (np) interaction remains unchanged at linear order in 5 np . In 
general, we expect these corrections to be small, though the difference between the in-medium 
pp and nn interactions can be important as we discuss in the next section. 

IV. NUMERICAL RESULTS 
A. Partial wave projection 



For the purpose of demonstration and for nuclear structure calculations it is convenient to 
evaluate the above in-medium NN interaction in the partial-wave basis consisting of eigenstates 
with definite L, S, J quantum numbers. We have followed the detailed description in ref. [43 ] 
and obtain with Uk = Vk + (4/ — 3)Wk , (K = C, S, T, SO, Q) the relevant linear combination 
in a state with total isospin / = 0, 1 the following projection formulas: 
a) Singlet matrix element with S = and L = J: 



(J0J\V NN \J0J) = - I (/: 
2 7_i 



b) Triplet matrix element with S 

(J1J\V NN \J1J) = ; 



U c - SU S 
1 and L = J: 



q 2 U T + p\z 2 - 1)U Q 



P.i(z) 



-i 



dz\ 2p z 



U so -U T + P 2 zU q 



{Pj +1 (z) + Pj^(z) 



+ 



U c + U s + 2p 2 (l + z)U T - Ap 2 zU so - p 4 (Sz 2 + l)U Q 



PAz) 



(3? 



(39) 
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c) Triplet matrix elements with S = 1 and L = J ± 1: 



(J±1,1J\V NN \J±1,1J) = - I dz{2p 2 



U SO ± 



2J+1 



U T -p 2 zU t 



Q 



PAz) 



+ 



U c + U s + p 2 [p 2 (l - z 2 )U Q - 2zU so ± 



P 2 U Q - U T 



2J+1 

d) Triplet mixing matrix element with S — 1, L' = J — 1 and L = J + 1: 

\fJTTp 2 



(40) 



[J-1,1J\V NN \J+1,1J) 



+ 



' 2 -J^dz{(u T -p 2 U Q )Pj +1 (z) 
(2 J - z{2 J + l))U T + P 2 zUq\ Pj(z)} . 



V1(2J + 



(41) 



Here, Pj(z) are ordinary Legendre polynomials of degree J and the momentum transfer q is set 



to q = py 2(1 — z). The resulting partial wave amplitudes are now functions of the momentum 
p and kf (or equivalently the nucleon density p). Note that the off-diagonal mixing matrix 
elements arise exclusively from the tensor operator d\ ■ qo<i ■ q and the quadratic spin-orbit 
operator B\ • (q x p) a 2 ■ (q x p). 



B. 5- waves and S — D mixing 

In Figs. H] and [5] we show by the solid lines the (diagonal) partial-wave amplitudes of Vj ow _k 
in the channels 1 So, 3 Si, and 3 Di as well as the 3 Si— 3 Di mixing matrix element for a cutoff 
Aiow-k = 2.1 fm _1 . We also show the modification of the low- momentum interaction due to the 
six different components of the density- dependent NN interaction derived from the leading-order 
chiral three- nucleon force in an isospin symmetric nuclear medium with density p = p$. The 
pion self-energy correction V^ 1 ' 1 and the Pauli-blocked one-pion exchange vertex correction 
V^r d ' 2 , which both arise from the long-range three- nucleon force, vanish as q — > and therefore 
give no modification to the S- wave amplitudes as p —>■ 0. Moreover, to a large extent they 
contribute with opposite sign and equal strength. Consequently, their combined effect on the 
^-wave amplitudes is rather small. In comparison, the long-range Pauli-blocked 27r-exchange 
term V^ d ' 3 contributes strongly to the S waves at large distances (small momenta) and has 
little effect at short distances (large momenta). In fact, Pauli-blocked 27r-exchange is the 
dominant mechanism for suppressing the S"-wave attraction at low momenta. From Fig. [5] 
we observe that the contributions V^ d ' 4 and V^f' 5 , resulting from the medium-range three- 
nucleon force proportional to cp, are naturally small in S- waves. Additionally, they enter with 
opposite sign and therefore play a very minor role in modifying the S'-wave interaction. The 
final contribution V^ 1,6 comes from the three- nucleon contact interaction proportional to eg. 
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FIG. 4: Modifications to the 1 So, 3 Si, 3 Di, and 3 Si— 3 Di partial wave amplitudes of V| ow -k (shown 
by the solid curve) due to the first three density-dependent contributions V™^' 1,2,3 at nuclear matter 
saturation density po = 0.16 fm -3 . 

Since it is independent of the momentum p (aside from the suppression above p = 1.5 fm -1 
due to the regulator function), it affects only the two S- waves. In combination with V^ d ' 3 it 
plays a major role in reducing the S-wave NN attraction in the nuclear medium. 

In the 3 D! partial wave we find again that the pion self-energy correction V™^' 1 largely can- 
cels the effect from the long-range Pauli-blocked vertex correction V^ d ' 2 . The Pauli-blocked 
27r-exchange contribution V^ d ' 3 now vanishes at small momenta and has a large effect at high 
momenta, in contrast to its role in the S-waves. By itself, V^r d ' 3 would nearly double the 
repulsion in the 3 Di channel. This enhancement of the 3 Di repulsion is reduced in part by the 
attractive contribution coming from the medium-range Pauli-blocked one-pion-exchange vertex 
correction V^ d ' 4 . The final possible contribution, V^ d ' 5 , in this channel is neglegible. Finally, 
the 3 Si— 3 Dx mixing matrix element receives contributions from all density-dependent terms 
except the short-range contact interaction V^r d ' 6 since each of the first five contributions in- 
cludes a tensor component. Only V^ d ' 2 and V^r d ' 4 give repulsive contributions to the 3 Si— 3 Di 
mixing, and these terms are generally weaker than the three attractive contributions. 

In Fig. [6] we show the modifications of the momentum-space matrix elements in the 1 So, 
3 Si, and 3 D X partial waves of V[ ow _k due to all density- dependent contributions. The complete 
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FIG. 5: Modifications to the 1 So, 3 Si, 3 Di, and 3 Si— 3 Di partial wave amplitudes of Vjow-k (shown 
by the solid curve) due to the three density-dependent contributions V^ N ' ' 5 ' 6 at nuclear matter 
saturation density po = 0.16 fm~ 3 . 

interaction is shown at the densities po/2 and po- The overall effect is to decrease the strong 
attraction in S- waves and to increase the repulsion in the (diagonal) 3 Di channel. The 3 Si— 3 Di 
mixing matrix element becomes on average only mildly more attractive. The repulsive effects 
increase with the nuclear density and this way provide the mechanism for nuclear matter 
saturation when using the low-momentum interaction V[ ow _k [h]] • Finally, in Fig. [7] we plot 
the in-medium pp interaction at nuclear matter saturation density po for isospin asymmetries 



& 



rip 



±0.2, where 5 np = 0.2 corresponds to the asymmetry reached in heavy nuclei. For 



larger isospin asymmetries, the effects scale linearly with S np . We plot as well the results for 



'rip 



—0.2, which corresponds to the nn interaction in a medium with 5 np = 0.2. As discussed 



in Section IIIIBI the effects on the pp interaction and nn interaction are equal and of opposite 
sign. Although the effects due to an isospin asymmetry are in general small, the difference 
between the pp and nn interactions can be significant, particularly in the 3 Si channel. 
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FIG. 6: Dependence of the low- momentum in- medium NN interaction on the nuclear density p. Shown 
are the momentum space matrix elements in the 1 So, 3 Si, and 3 Di partial waves and the 3 Si- 3 Di 
mixing matrix element. 



P [fm"'] p [fm -1 ] 
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FIG. 7: Dependence of the low-momentum in- medium nucleon-nucleon interaction on the isospin 
asymmetry S np = (p n — p p )/p at po = 0.16 fm~ 3 . Shown are the momentum space matrix elements in 
the 1 So, Si, and 3 Di partial waves and the 3 Si- 3 Di mixing matrix element. 
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FIG. 8: Modifications to the 1 Pi, 3 Po, 3 Pi, and 3 P2 partial wave amplitudes of Vi ow -k (denoted by 
the solid line) due to the first three density-dependent contributions ' 1,2,3 at a nuclear density 
P = Po- 

C. P- waves 

There are three uncoupled L = 1 partial waves, namely 1 Pi, 3 Po, and 3 Pi, as well as the 
3 P2 partial wave which can couple through the tensor and quadratic spin-orbit forces to the 
3 F2 partial wave. In Figs. [8] and [9] we have plotted the effects of the six components of the 
density-dependent NN interaction on each of these P-wave amplitudes. We consider a medium 
of isospin-symmetric nuclear matter at saturation density po- Again we find that although 
the pion self-energy correction and the long-range one-pion-exchange vertex correction are the 
largest contributions in all channels, taken together they have only a moderate effect. The Pauli- 
blocked two-pion exchange contribution V^ d ' 3 is attractive in the P-wave spin-singlet states 
and repulsive in spin-triplet states, though only in the 3 Pi and 3 P2 does it play a significant 
role. Effects from the one-pion exchange vertex correction arising from the medium-range three- 
nucleon force (V^ d ' 4 ) can be as large as 20%, but those from the vertex-corrected short-range 
nuclear force (V^f' 5 ) are negligible. 

In Fig. [in] we have plotted the complete P-wave interactions at densities po/2 and po- Aside 
from the 3 Pi channel, which on average receives a small attractive contribution from the sum of 
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FIG. 9: Modifications to the Pi, Po, Pi, and P2 partial wave amplitudes of Vi ow _k (shown by the 
solid line) due to Vjfffl ' 5 at a nuclear density p = po. 



the six density-dependent terms, we see that the remaining L = 1 partial waves all receive a net 
repulsive contribution. The 3 Po partial wave is particularly sensitive to these modifications; at 
nuclear matter saturation density, nearly all of the attraction at small momenta vanishes and 
the repulsion at higher momenta is largely increased. This results primarily from the repulsion 
due to V^ d ' 2 , V^ d ' 3 , and V^ d ' 4 . The 3 Pi and 3 P2 partial waves are less sensitive to the density 
dependent terms. The 3 Pi channel receives two attractive (V^ d ' 2 , V^f' 4 ) and two repulsive 
(^jvjv' 1 ; ^JVjv d ' 3 ) modifications which give rise to only a small net repulsive effect limited to 
intermediate momenta. Finally, for the 3 P 2 partial wave, only V^f' 1 , V^ d ' 2 , and V r ™^ d ' 3 are 



3 1 

rmed,2 



NN 



important. The net repulsion from V™^ z and ' 3 is approximately twice as large as the 
small attraction from V™^' 1 - Isospin asymmetry effects (5 np = ±0.2) in the L = 1 partial waves 
are shown in Fig. [TU The largest effect is in the 3 Po channel where the difference between the 
pp and nn interaction is approximately 25% of the total interaction strength. In fact, in the pp 
channel there is almost no attraction remaining. 
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FIG. 10: Dependence of the low-momentum in-medium nucleon-nucleon interaction on the nuclear 
density p. Shown are the momentum space matrix elements in the Pi, 3 Po, 3 Pi, and 3 P2 partial 
waves. 
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FIG. 11: Dependence of the low-momentum in-medium nucleon-nucleon interaction on the isospin 
asymmetry 5 np = (p n — p p )/p at po = 0.16 fm~ 3 . Shown are the momentum space matrix elements in 
the ^i, 3 P , 3 Pi, and 3 P 2 partial waves. 
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D. D- waves and D — G mixing 



Although the remaining L = 2 partial waves are mostly smaller by an order of magnitude 
than the L = and 1 partial waves, it is important to see whether the trends observed in 
the lower partial waves continue. In Figs. [12] and [13] we have plotted the effects from the five 
(nonvanishing) density- dependent contributions in the 1 D 2 , 3 D 2 , and 3 D 3 partial waves. As in 
the S and P-waves, the V^f' 1 and V^ d ' 2 terms contribute strongest and with opposite sign. 
The long-range Pauli-blocked 27r-exchange contribution V^ d ' 3 has little effect, except in the 
X D 2 partial wave where it gives rise to a strong repulsion. Together with V^ d ' 2 , the Pauli- 
blocked 27r-exchange contribution gives rise to a repulsion in this channel that is significantly 
larger than the attractive contribution from the pion self energy term V™^' 1 . In fact, both 
V^f' 2 and V^r d ' 4 are repulsive in all of these channels, except in the 3 D3 partial wave where 
they combine to generate an overall attraction. These results can be seen from Fig. [TH in which 
we have plotted the complete interaction at densities po/2 and po- From this figure one can 
see that all L = 2 partial waves, except the 3 D 3 channel, become less attractive due to the 
density-dependent corrections. In Fig. [15] we show the effect of an isospin asymmetry on these 
partial wave amplitudes for 5 np = ±0.2. Again, we find that in general the modifications are 
small, though in some channels, such as the difference between the pp and nn interactions 
cannot be neglected. 
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FIG. 12: Modifications to the 1 D2, 3 D2 and 3 D3 partial wave amplitudes and the 3 D3— 3 G3 mixing ma- 
trix element of Mow-k (shown by the solid line) due to the first three density-dependent contributions 
1,2,3 at saturation density po- 



P [fm" 1 ] p [fm" 1 ] 




P [fm [ ] p [fm 1 ] 

FIG. 13: Modifications to the 1 D 2 , 3 D 2 and 3 D 3 partial wave amplitudes and the 3 D3— 3 G3 mixing 
matrix element of Vi ow -k (shown by the solid line) due to V^ d ' 4 ' 5 at saturation density po- 
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FIG. 14: Dependence of the low-momentum in-medium nucleon-nucleon interaction on the nuclear 
density p. Shown are the momentum space matrix elements in the 1 D2, 3 D2 and 3 D3 partial waves as 
well as the 3 D 3 - 3 G 3 mixing matrix element. 
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FIG. 15: Dependence of the low-momentum in-medium NN interaction on the isospin asymmetry 
3np = (Pn — Pp)/p & t Po = 0.16 fm -3 . Shown are the momentum space matrix elements in the 1 T>2, 
3 D2 and 3 D 3 channels and the 3 D 3 — 3 G 3 mixing matrix element. 
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V. SUMMARY AND CONCLUSIONS 



Using chiral effective field theory, we have derived in this work density-dependent corrections 
to the nucleon-nucleon interaction in isospin-symmetric as well as isospin-asymmetric nuclear 
matter. These corrections have been calculated from the six one-loop in-medium NN-scattering 
diagrams generated by the leading-order chiral three-nucleon interaction. The resulting in- 
medium NN interaction has been transformed into the \LSJ) basis. Although we have combined 
the density-dependent corrections with the low-momentum NN potential Viow-k at a selected 
cutoff Ai ow _k = 2.1 fm _1 , the analytic expressions for the density-dependent terms should hold 
at any resolution scale by re-adjusting the values of the running low-energy constants co(A low _ k ) 
and cs(Ai ow _k). In the present work the values of the two low-energy constants cd = —2.06 
and ce = —0.63 were taken from a previous calculation 38] of A = 3,4 binding energies using 
the low-momentum NN potential Vi ow -k at Ai ow _t = 2.1 fm _1 . 

After projecting the interaction into partial waves, we find that the largest density-dependent 
corrections come from the long-range Pauli-blocked pion self-energy and vertex correction. 
However, to a large extent they cancel in all partial waves studied and thus give rise to effects 
comparable to the remaining density-dependent corrections. The two contributions from the 
medium-range three-nucleon force proportional to cd are small and generally play only a minor 
role. In contrast to this, the long-range Pauli-blocked 27r-exchange as well as the contact inter- 
action proportional to ce provide significant repulsion in most partial waves. The latter acts 
only in 5-waves where it decreases the attraction in the spin-singlet and spin-triplet channels 
by approximately 20%. In fact, the net effect of all density-dependent corrections is repulsive 
in nearly all partial waves channels we studied. The repulsive effects increase with density and 
in this way provide a mechanism for nuclear matter saturation [h]]. The repulsive (stabilizing) 



nature of the pion-induced three-body 



calculation of nuclear matter in ref. 



3l| 



orces is also observed in the chiral perturbation theory 



Furthermore, we have found that additional medium modifications due to a (small) isospin 
asymmetry are present only in the pp and nn channels and are generally small. However, the 
resulting difference between the in-medium pp and nn interactions can be significant in certain 
partial waves. This work should serve as a foundation for future nuclear structure calculations 
exploiting density-dependent chiral effective two-nucleon interactions. 
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